Length-Dependent Conductance of a Spin-Incoherent Hubbard Chain: Monte Carlo 

Calculations 



o 
o 

(N 
5h . 

< 
00 



Olav F. Syljuasen 1 
1 NORDITA, Blegdamsvej 17, DK-2100 Copenhagen 0, Denmark 

The dc conductance of a short spin-incoherent Hubbard chain coupled to leads is investigated 
using quantum Monte Carlo calculations. In contrast with the Luttinger liquid regime, where the 
conductance is equal to the non-interacting value, the spin-incoherent regime displays a conductance 
that decreases rapidly with chain length down to a value of roughly 1.5e 2 /h for a four site chain 
followed by a slower decrease for longer chains. We also discuss the resistance contribution from 
scattering in the contacts. 

PACS numbers: 71.10.Pm, 73.23.-b, 71.10.Fd, 73.63.Nm 
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Advances in making small electronic devices have made 
it increasingly important to understand the quantum 
mechanics of confined electrons. Especially important 
is the quantum wire where the electrons are restricted 
to move in a narrow channel. When such a quantum 
wire is made thinner its conductance decreases in per- 
fect steps of 2e 2 //i0]. There is, however, an anomaly 
in this otherwise perfect picture; at 0.7 x 2e 2 /h an ex- 
tra steplike structure is seen at very low electron den- 
sities at finite temperatures 0. While the regular steps 
can be explained in terms of non-interacting electrons, 
different theories have been proposed for the anomalous 
stepd |, i, S]. A particularly intriguing theory[3| is 
that the suppressed conductance is a result of the elec- 
tron gas entering a spin-incoherent (SI) regimeQ at high 
enough temperatures. It is however unclear how the pref- 
actor 0.7 comes about within this theory. Experiments 
indicate that this prefactor depends on the length of the 
wire[8|. However, not much is known about this theoret- 
ically. The aim of this Letter is to investigate this length 
dependence numerically. 

At zero temperature, T = 0, the low-energy model 
of a generic interacting ID electron liquid is the Lut- 
tinger liquid (LL) which is characterized in part by hav- 
ing independent spin and charge excitations each hav- 
ing separate coupling constants, characteristic velocities 
and bandwidths. The conductance of a uniform pure 
LL depends on the charge sector coupling constant K c 
alone and is at low temperatures G = 2K c e 2 jh% [ijj. 
However this result does not coincide with what is being 
measured in experiments where leads are inevitably con- 
nected to the interacting electron liquid. The leads have 
a profound influence on the dc conductance and it has 
been shown by modeling the leads as ID non-interacting 
electron gases that the conductance remains at the non- 
interacting value 2e 2 /h provided one couples the leads 
sufficiently smoothly to the interacting wire fill]. Other- 
wise the conductance is solely determined by scattering 
in the contacts. As will be shown the presence of leads 
also influences the conductance in the SI regime. 

For temperatures higher than the bandwidth of ei- 



ther the spin or the charge excitations one do not ex- 
pect the LL description to hold. In general the spin 
bandwidth J is smaller than the charge bandwidth Et 
when the interaction energy dominates. Thus a situation 
where J <C T <C Ef is conceivable. This is known as 
the SI regime, where in contrast to the LL regime, the 
electron Green function displays non-propagating spin 
excitations 0, E2], a broad momentum distribution [3], 
and an anomalous density of states 14 , li| . The bound- 
ary Green function has also been explored (ljj as well as 
Coulomb drag effects 17j , the Fermi edge singularity UM 
and transport properties in the presence of impurities [l9(. 
It was argued in Ref. |3|] that spin-charge separation gets 
violated when leads are coupled to a SI wire, and as a con- 
sequence the dc conductance is renormalized to a value 
e 2 //iQ. While the result in Ref. 0] was obtained for an 
electron liquid at low densities forming a Wigner crys- 
tal, it has been argued that the physics of the SI regime 
is largely independent of microscopic details, and that 
one might as well replace the electron gas with a model 
having short-range interactions like the Hubbard model. 
This was utilized in Ref. [201 ] where the non-equilibrium 
conductance and noise of a SI chain coupled to leads was 
calculated in two regimes, a strongly biased regime and 
a regime with dissipative spin damping, both giving a 
conductance e 2 /h. 

We model the SI wire coupled to leads as a Hubbard 
model with a site dependent potential Ui 



(ij>,<r={Uj 



C ia C j<? 
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where ci a is the spin a fermion annihilation operator at 
site i, and is the corresponding density operator. 
We treat a ID lattice with L sites using open bound- 
ary conditions and divide it into five regions, see Fig. [1] 
Two non-interacting lead regions each of length Ll where 
u = 0, an interacting chain region where u — U and 
two contact regions each of length Lq where u is site- 
dependent so that the particle density in the chain region 
interpolates roughly linearly to the density in the lead. 
Throughout this Letter we add a uniform chemical po- 
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FIG. 1: The ID lattice with noninteracting lead regions of 
length Ll, contact regions of length Lc, and the fully inter- 
acting chain of length Lw ■ 

tential fx = —0.3t causing the density in the lead regions 
to be ~ 0.9, slightly less than half-filling. The tempera- 
ture is T = t/80 and L = 192. To avoid complications 
due to Kondo physics, we keep L\y even, yet we expect 
that our results also apply to odd Lw above the Kondo 
temperature. 

For the Hubbard model the bandwidth of the spin ex- 
citation spectrum is J = (At 2 /U)n[l — sin(27rn)/27rn], 
while the bandwidth of the charge spectrum is propor- 
tional to the Fermi energy Ef. The SI regime J <C Ef 
is thus realized for low densities, or equivalently for large 
U . In addition T must be placed in-between these scales. 
The latter is more demanding for low densities than for 
large U, because of the smallness of Ef at low densities, 
thus requiring very low temperatures which are costly 
using the numerical technique employed here. Therefore 
we will model a chain in the SI regime by setting U = oo 
in the interacting region. This implies J = even for 
high densities, thus Ef can be kept of the order t. 

In order to investigate the ID inhomogeneous system 
Eq. (Q} we employ the Stochastic Series Expansion (SSE) 
quantum Monte Carlo (QMC) method with directed-loop 
updates [lH. Fermion QMC simulations frequently come 
with a minus-sign problem. However this is avoided in 
ID using open boundary conditions. The fermions are 
represented in the occupation number formalism follow- 
ing the sign convention described in Ref. 12211 , and the 
directed-loop rules were taken from Ref. [23j. In addi- 
tion to having updates that add/remove a particle with 
a certain spin we also employ moves that flips the spin 
of a particle. This is necessary to ensure short autocor- 
relation times in the SI regime. 

The dc conductance G is evaluated as the 
linear response of the current operator j x = 
i{et/fi)J2a=^i( c ia c x+ia - 4+i<r c ^cr) at position x 
to a discontinuity in the chemical potential at y in the 
limit of zero frequency 

i r°° 

G=lim 9 (z), g { z )=Re- dte izt (\j x (t),P y }) (2) 

where P v is the sum of fermion charge density opera- 
tors at sites to the right of y, P y — eJ2 y >> y n y'- The 
extrapolation to zero frequency can be taken along any 
path in the complex plane. We will extrapolate along the 
imaginary axis[24J, thus z is taken to be imaginary and 
is denoted z = ito. Using the charge/current continuity 
relation for a ID system with open boundary conditions 
one finds that g(iu>) can be evaluated at the Matsubara 
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FIG. 2: (Color online) Imaginary frequency conductance for 
free fermions. The dotted, dashed and long dashed lines are 
exact results for different system sizes L indicated by the leg- 
ends. The circles are QMC data at the three lowest non-zero 
Matsubara frequencies for L = 192. The solid lines are the 
extrapolated QMC results using the prescription in the text 
and the error bars at u = indicate the dc conductance in 
the thermodynamic limit. The upper set of curves (black) is 
for x = y = L/2 — 1 and the lower (blue) is for x = L/2 — 1, 
y — x + 1. The three vertical bars indicate the frequency 
uj = 2nhvf / L for the three system sizes L = 768, 384, 192. 

frequencies uj = lu„ = 2imT as[25j 

g(i(J n ) = -^j d,TCOs(0J n T){P x (T)P y (0)}. (3) 

This correlation function which involves only density op- 
erators is easily evaluated using the QMC method. While 
for a finite system, g(iuj = 0) = 0, the correct way of ob- 
taining the conductance in the thermodynamic limit is to 
first take the infinite system size limit and then extrapo- 
late to zero frequency. For a big enough system the finite 
Matsubara frequencies is not appreciably affected by the 
system size, thus the dc conductance in the thermody- 
namic limit can be gotten from extrapolating the conduc- 
tance at the finite Matsubara frequencies to zero along 
the imaginary frequency axis. Figure [3 shows results for 
free fermions. The two set of curves are for different 
choices of x and y. Both lead to a value consistent with 
the exact value 2e 2 /h. The circles are QMC data and 
the extrapolations are gotten by constructing a rational 
polynomial function of degree [p/q] that coincides with 
the QMC data at the lowest p + q+1 Matsubara frequen- 
cies. As the QMC data has error bars this construction is 
bootstrapped over 10 4 repetitions, and the median value 
is recorded. The median is chosen instead of the average 
to minimize the effects of spurious poles. This is repeated 
for all [p/q] with 5 < p + q < 8 and p, q > 2. The solid 
line shows the average of these (median) values and the 
error bar at lu — shows the corresponding maximum 
spread, which is larger than the error of a single [p/q] 
extrapolation. 

The dashed and dotted curves that approach for low 
lu are exact results on chains of different lengths. They 
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FIG. 3: (Color online) Imaginary frequency conductances ex- 
trapolated to zero frequency for uniform Hubbard chains with 
different values of U. From top to bottom U/t = 0, 1, 2, 4, oo. 
x = L/2 — 1. The circles are for y = x and the squares are 
for y = x + 1. 



start to deviate from the extrapolated results at u « 
hvf2ir/L which correspond to frequencies lower than the 
level spacing of the system. This sets a lower limit, L > 
hvf/3, for the system size needed at a finite temperature 
T = 1/(3, in order to ensure that the finite Matsubara 
frequencies all reflect the infinite size limit. For T = t/80 
a system size of L = 192 is sufficient. 

To check our method we obtain the dc conduc- 
tances for uniform chains with no leads having differ- 
ent values of the interaction U. The measured parti- 
cle densities are 0.903, 0.781, 0.696, 0.601, 0.457 for U/t = 
0,l,2,4,oo respectively. Fig. [3] shows the results con- 
taining two different extrapolations for each value of U. 
Taking the average of the extrapolated results we ob- 
tain G/(e 2 /h) = 1.999 ± 0.004, 1.817 ± 0.008, 1.656 ± 
0.012, 1.431 ± 0.006, 0.999 ± 0.014 for U/t = 0, 1, 2, 4, oo 
respectively. These values agree well with G = 2K c e 2 /h, 
where K c is the Luttinger liquid charge sector coupling 
constant whose value can be obtained for the Hubbard 
model for a given density and U by the Bethe Ansatz[2fj|. 

We now attach leads to the U = 4 Hubbard chain. 
Figure |4] upper panel shows the imaginary frequency con- 
ductances and their extrapolations to zero frequency for 
different lengths Lw of the interacting region. In order 
to make the contact resistance small we have made con- 
tact regions of length Lc = 2. While at high Matsubara 
frequencies the conductance curves are close to the con- 
ductance curve of the uniform chain without leads, there 
is a change in behavior when the frequency gets below 
lj ~ hvf/L\Y, where the conductance curves extrapo- 
late towards the non-interacting value in accordance with 
Refs. 
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To see how the value of the dc conductance 
depends on the length of the contact regions we show in 
Fig. |4] lower panel the conductance for different lengths 
of the contact regions keeping the length of the interact- 
ing chain fixed. The curve with an abrupt contact region 
(Lc = 0) shows a lower dc conductance than the others 
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FIG. 4: (Color online) Imaginary frequency conductances for 
interacting (U = At) chains with leads. Upper panel: The 
curves are for different Lw indicated by the legends and a 
fixed Lc = 2. For comparison the uniform U — At chain 
without leads is also shown (labeled oo). Lower panel: The 
curves are for different Lc indicated by the legends and a 
fixed Lw = 8. In both panels x — y = L/2 — 1. 



which appear to have reached the adiabatic limit already 
for Lc = 1. 

The SI limit, U — oo, is reached by restricting the 
Hilbert space so that no sites in the interacting region 
are doubly occupied. Fig. [5ja) shows the imaginary 
frequency conductances and their extrapolations for SI 
chains of different lengths Lw coupled to leads. As in 
the LL regime, the data at low imaginary frequency differ 
from that of the uniform chain without leads (labeled oo) 
and extrapolate to values larger than e 2 /h. The extrap- 
olation for the longest chains are omitted because only 
the few lowest Matsubara frequencies distinguish these 
curves from the uniform one causing excessively large er- 
ror bars. Nevertheless it is quite clear that the dc con- 
ductances of the longest chains also extrapolate to values 
larger than e 2 /h. This is rather remarkable in view of the 
exponential decrease of the single-particle Green function 
in the SI regime. The extrapolated dc conductances are 
plotted in the inset, one datapoint for each \p/q], and 
reveal a rapid decrease with increasing chain length to a 
value ~ 1.5e 2 //i for Lw — 4, and then a further slower 
decrease. Note that this intermediate conductance value 
is close to the value relevant to explain the 0.7 conduction 
anomaly. 

In addition to the intrinsic resistance of the interact- 
ing region, our results include possible contributions from 
scattering in the contact regions. The contact regions 
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FIG. 5: (Color online) Imaginary frequency conductances of 
SI chains with leads, y — x — L/2— 1, L — 192. (a) curves are 
for different Lw indicated by the legends, Lc = 2. The inset 
shows the dc conductances, (b) Fixed Lw = 8 and different 
Lc indicated by legends. Inset: dc conductances, (c) particle 
density vs. site for different Lc- Lw = 8. The Lc ~ 4(16) 
curve has been shifted downwards by 0.1(0.2) for clarity. 



used were of length Lc = 2, which is long enough to 
give adiabatic contacts in the LL regime, see Fig. 2J To 
check their role in the SI regime we show in Fig. |5£b) 
results for different Lc- While there is a significant re- 
sistance contribution for abrupt junctions, Lc = 0, there 
are only small differences between the dc conductances 
for contacts of lengths Lc > 1, see inset. It is not entirely 
clear however, due to the large error bars for the largest 
values of Lc, if the dc conductance has saturated or will 
keep increasing for even larger Lp's. Assuming satura- 
tion we estimate that the contact contribution to the re- 
duced conductance is ~ 0.1e 2 /h for Lc — 2. The density 
variations in the contacts causing scattering is shown in 
Fig. [5fc) . Note the rather abrupt density changes that 
cause resistance in the case Lc — 0, compared to the 
smoother density variations for Lc = 16. 

Our results show that the dc conductance decreases 



then decreases further slowly with Lw for Lw > 4. It is 
not clear what sets the rather long length scale associated 
with this slow decrease. It is plausible that this is related 
to the magnitude of terms in the effective Hamiltonian 
breaking spin-charge separation. While it is intriguing 
that the conductance values obtained here are consistent 
with the value 0.7 x 2 for a broad range of wire lengths 
Lw, it remains to be seen to which extent these values 
are universal. 



The simulations were in part carried out using com- 
puters provided by the University of Aalborg, Denmark 
and the University of Tartu, Estonia using the NorduGrid 
ARC middleware. 
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